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Abstract. The problem of calculating the period of second order nonlinear 
autonomous oscillators is formulated as an eigenvalue problem. We show that 
the period can be obtained from two integral variational principles dual to 
each other. Upper and lower bounds on the period can be obtained to any 
desired degree of accuracy. The results are illustrated by an application to the 
Duffing equation. 



1. Introduction 

The motion of a particle in one dimension under the action of a conservative force 
is described by the simple equation u = —f{u), where the force — / depends on the 
amplitude of the motion u and one or more parameters. For such conservative sys- 
tems the total energy E = u 2 /2 + V(u), is conserved. While a complete description 
of the motion can be obtained from the analysis of phase space, an explicit analytic 
expression for the period of oscillations around equilibria is not always available. 
For linear problems the period does not depend on the amplitude of the motion, 
but when nonlinearities are present, even though we know the expression for the 
period, namely 

du 



(1) T = 2 



y/2(E-V(u)) 



where uq and u\ are turning points, the integral cannot be evaluated analytically 
except for special nonlinear potentials V{u). The usual procedure in such cases is 
to consider the nonlinearities as small perturbations to the linear problem and to 
obtain a series expansion for the period of oscillation. To obtain a close estimate 
of the period at large amplitude, the perturbation expansion must be carried to 
high orders. Often the series may be divergent. Recently a method to obtain a 
rapidly convergent series has been developed 0, based on Kleinert's variational 
perturbation theory (see Chapter 5 of 0). 

In previous work we studied a nonlinear eigenvalue problem, of which the above 
problem constitutes a particular case, and showed that the eigenvalue (the period in 
the present case) can be obtained from an integral variational principle. From it, we 
obtained arbitrarily accurate upper bounds on the period, valid at all amplitudes, 
without resource to perturbation theory. The purpose of the present work is to show 
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the existence of a dual variational principle from which lower bounds arc obtained, 
thus providing a complete variational characterization of the period. 

We consider systems with an equilibrium point u = 0, and will study the period 
of oscillations around this equilibrium. Moreover, for the sake of simplicity, we 
assume that the force is an odd function of u. The results can be generalized in a 
simple way to a general force term. If the force is odd, the period can be evaluated 
by considering the motion through a quarter of a period. Chosing the quarter 
period in the quadrant {u < 0, u > 0) of phase space, the period may be calculated 
solving il — —/(it), with u(0) = u m , ii(0) = and it(T/4) = 0. 

Introducing the new variable r = At/T, the problem reduces to finding the 
eigenvalue A of 

(2) -u" = \f(u), u(0) = u m , it'(0) = 0, u{l) = 0. 

The period is given by T = 4\/A. Here, primes denote derivatives with respect to 
t. Notice that the eigenvalue A depends on the amplitude u m - 

In previous work we proved that the eigenvalue A is characterized by the varia- 
tional principle 

(3) A[u m ] = max — u — , 

9 2 Jo m f(u)g(u)du 

where the maximum is taken over all positive functions g(u) such that g(0) = 0, 
g'(u) > 0. The maximum is achieved for g — g which satisfies 

9W {E-V{u)f/f 
An application of this result to the Duffing equation yielded a simple approximate 
formula as a close upper bound for the frequency uj — 2tt/T. In the following we 
construct the dual of The two variational principles give a complete charac- 
terization of the eigenvalue problem. In Section 2 we give the complete variational 
characterization of the period. In Section 3 we apply the variational principle and 
its dual to estimate the period of the Duffing oscillator. 

2. A VARIATIONAL PRINCIPLE AND ITS DUAL 

Many authors have considered the following two point boundary value problem 

(4) -y" = Xf(y) in (0, 1) 

with 

(5) 2/(0) = 0, 2/(1) = 0, and y(0) = y m . 

Here f(y) is positive and continuous but not necessarily when y = 0. For example, 
this type of nonlinear two point boundary value problem arises in the study of heat 
generation and stability of temperature distribution of conducting plates |SJ . Under 
the assumption that / is positive and continuous, Laetsch [5] proved the existence 
and regularity of nontrivial solutions to In fact, the positive solutions to 

©,© are decreasing and the problem can be reduced to a quadrature, 

(6) \y' 2 + V(y) = E, 
where V(y) = A f(s) ds and E = V(y m ). 
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In 0], Theorem 2.7, we proved the following variational characterization for the 
principal eigenvalue of ©,0. 



Theorem 1. Let the pair (X,y) be the principal solution (i.e., with y(x) > 0) of 
the two point boundary value problem 

dx 2 

subject to y'(0) = y(l) = 0. Let y m = y(0) be the sup-norm of the solution. Here, 
f(y) is a general nonlinear term which is positive and continuous in (0,y m ). Then, 



(7) -33 = A/(») 



(8) A [ym] = max - ' , 

geD 2 Jo f(y)9{y)dy 

where D = {g | g G C 1 (0, y m ),g' > 0,5(0) = 0}. Moreover, the maximum is 
attained at a unique (up to a multiplicative constant) g G D. The maximizing g 
satisfies 

(9 > I - (*-*(,»"■ °" ,i! ' (0) = - 

In what follows it is convenient to introduce 

(10) h{y) = F(y m ) - F(y), 
where 

(11) F(y)= I" f(s)ds. 



Remark 1. Using the optimal g, given through equation 0), in integrating the 
denominator by parts, one gets the value of X[y m ] for the principal solution. This 
is given exactly by 

,i2 » *"HCTtM 2 ' 

After small manipulations, this expresion for X[y m ] is precisely the formula for the 
peiod given by 0) in the introduction. 

In the sequel, we will obtain a new variational characterization of A[y m ] which 
is the dual of JSJ . The dual principle is given by the following theorem. 

Theorem 2. Let the pair (X,y) be the principal solution (i.e., with y{x) > 0) of 
the two point boundary value problem Then, 



1 2 J*»(l -2(r(»))VW») + l)/h(s)ds 

— = max 
A <jeh 



(13) — = max 



s 

where H = {a | a G C(0, y m ), & > 0}. The maximum is attained at a = and it is 
given precisely by 

,i4) j-{r7w/Y- 
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Proof. We will start from the original variational characterization of A embodied 
in Theorem 1, and use the standard Fenchel-Moreau duality (see, e.g., the book 
by Ekeland and Temam ) to obtain our dual principle. In order to simplify our 
calculations it is convenient to define q(y) — g'iy) 1 ^ 3 , with g 6 D. Integrating the 
denominator of JSJ by parts, we obtain 

1 fn "* h(s)q 3 ds 
(15) -f-^ = 2min^ \Jl 

or, equivalently, 



3 ' 



1 



Vn 



(16) —. r = 2 min / h(s)q 3 ds, 

A[y m J qeD' J 

where D' = {q \ q G C{0,y m ),q > Q, J£ m q(s) ds = 1}. Let J[q] = f^ m q(s) 3 h(s)ds, 
and let J(t) be the Legendre transform of J, i.e., 

(17) J[t] = sup ( [ V t(s)q(s)ds - J[q] 

qeD' \Jo 

In principle we should compute J{t) for all t € L 2 (0,y m ). Doing so one finds that 
the function t that maximizes — Jit) is t = (we leave the details to the reader). 
Having this in mind, and for the sake of simplicity, we may restrict to nonnegative 
functions t(s), which we do in the sequel. In order to simplify the computation of 
the Legendre transform of J it is convenient to use the change of variables t — > a 
given by t = 3aer, with a given by 



(18) l = Va/ ^^iW *• 



tin 



' 1 + ojs) 
h(s) 

With this change of variables it follows from l|17|) and the definition of J[q] that 
(19) J\a] = «VS (2a- 




Using JTSJ in JTSJ finally gives, 

(2Q) _j w = ft(l-^(s))Vns) + l)/h(s)ds ^ 

Using Fenchcl Morcau duality, 

min J(q) = max —J (a) 

q a i 

and the theorem follows from here. 



□ 



3. An application: The Duffing equation 

To illustrate our variational principle and its dual, we will consider as an example 
the Duffing oscillator. Consider the two point boundary value problem 

(21) x= -X(x + 8x 3 ), 

on the interval (0, 1) with x(0) = and x(l) — 0. Denote by x m — x(0) the 
sup-norm of the solution. We are only interested in the principal branch (A, x m ), 



VARIATIONAL CALCULATION OF THE PERIOD OF NONLINEAR OSCILLATORS 5 



i.e., in the positive solution. In terms of A, the period of the oscillator is given by 
4yA. For this simple example one can compute the period of the oscillator in closed 
form in terms of complete elliptic integrals. However, since our purpose here is to 
illustrate the use of the dual variational principles, wc will use simple trial functions 
in Theorems 1 and 2 above to compute upper and lower bounds on the period of the 
Duffing oscillator. There is a vast literature on the Duffing oscillator, and the use 
of perturbative schemes to obtain the period of the oscilator (see e.g., [0] and the 
references therein). Most perturbative schemes yield divergent asymptotic series 
for the period which can, however, be converted into exponentially fast convergent 
series by Kleinert's variational perturbation theory |7|. 

To begin with, we use Theorem 1 to get a lower bound on A. As we pointed out 
in J3|, a good trial function is given by 

( 22 ) 9(x) = ± * = , 

which is certainly in D. This trial function will yield a lower bound on A which 
gives an excellent agreement with the exact value of A near the bifurcation point 
Ai = 7r/2. In fact, as it was illustrated in [I] it gives a good agreement for all values 
of the amplitude x m . From (|2"T|> we get 

(23) g'{x) ' 



04,-^)3/2- 

Using this trial function g in Theorem 1, yields the lower bound 
(24) A > ' 



4 1 + 3&4/4' 

If we denote by oj — 2tt/T, the above bound gives the following upper bound on to 
for the Duffing oscillator, 



(25) W <yi+~fe^. 

We now proceed to compute upper bounds on A (i.e., lower bounds on lu) using the 
dual principle given by Theorem 2. For the Duffing oscillator, f(x) = x + fa 3 , and 
therefore F(x) = x 2 /2 + fa 4 /4. Hence, we can write 

(26) h(x) = h {l + 6z) 
where 

(27) ho = ^(x 2 n - x 2 ), and z = ^{x 2 m + x 2 ). 

In order to get simple, computable, bounds using the variational principle given 
by Theorem 2, it is convenient to express the trial function a in terms of a new 
variable, p, through the following relation, 

(28) ?±1 = 1 >. 

n tiQ 

Different approximations (i.e., different upper bounds on A) are obtained by choos- 
ing different expressions for the trial function p. The simplest way to get computable 
bounds to different degree of accuracy is to use as trial p the Taylor expansion of 
(1 + fa) -1 / 2 truncated to any even power. In the sequel we give the details for 
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the first nontrivial truncation. Better approximations can be obtained in the same 
way. Choosing 

(29) p(z) = l-^+|5V, 

(which corresponds to the Taylor expansion of (1 + (5z) -1 / 2 truncated to second 
order), yields 

(30) a= ^f)!(40-15<5z + 9<SV), 
which is nonnegative, and therefore belongs to H. Denoting by 

A = 67108864 - 25165824^4 + 149422085 2 a4 - 41287680(5 3 x^ 

+31073280(5 4 x^ - 24748416<5 5 a; r 1 ° + 7192476(5 6 x^ - 2202957<5 7 x^, 

and 

B = 4(256 - 96^4 + 57(5 2 xf„) 3 , 

then, 
and 

(32) ^>\/J- 




0,0 0,2 0,4 0,6 0,8 1,0 1,2 1,4 



Figure 1. The solid line is the exact solution. The dot-dot-dash 
line is the upper bound, whereas the dashed lines are lower bounds 
from the variational principles. 
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In Figure 1, we plot the exact value (solid line) of u> = 2ir/T for the Duffing os- 
cilator as well as the upper bound (|25|l and the lower bound (|32|l . The upper bound 
is an excellent approximation for lj for all amplitudes, whereas more complicated 
trial functions have to be used to get better lower bounds onw. In the figure we 
also include the lower bounds obtained by truncating (1 + <5z) _1 / 2 to degree 4, 6, 8, 
and 10 as a function of z. It is interesting to point out that all these lower bounds 
are precisely the curves for lj obtained by standard perturbation theory. 
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